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Abstract 

We summarize a recent work on the subject title. The Dirac equation 
in a curved spacetime depends on a field of coefficients (essentially 
the Dirac matrices), for which a continuum of different choices are 
possible. We study the conditions under which a change of the coef- 
ficient fields leads to an equivalent Hamiltonian operator H, or to an 
equivalent energy operator E. In this paper, we focus on the standard 
version of the gravitational Dirac equation, but the non-uniqueness 
applies also to our alternative versions. We find that the changes 
which lead to an equivalent operator H, or respectively to an equiv- 
alent operator E, are determined by initial data, or respectively have 
to make some point-dependent antihermitian matrix vanish. Thus, 
the vast majority of the possible coefficient changes lead neither to an 
equivalent operator H, nor to an equivalent operator E, whence a lack 
of uniqueness. We show that even the Dirac energy spectrum is not 
unique. 

1 Introduction 

This paper summarizes a recent work [T] on the subject title. Thus, rather 
than to give detailed arguments and to present all relevant results, our aim 
here is to make the origin of this surprising result more easily visible. 
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Quantum effects in the classical gravitational field are observed, e.g. on 
neutrons [21 El H], which are spin | particles. This motivates work on the 
curved spacetime Dirac equation. The standard version of the latter (see 
e.g. Refs. [6] ) is due to Fock and to Weyl, and will be referred to as the 
Dirac-Fock- Weyl (DFW) equation. On the other hand, two alternative Dirac 
equations in a curved spacetime were derived recently [7J , by applying directly 
the "classical-quantum correspondence", i.e., the correspondence between a 
classical Hamiltonian and a quantum wave operator. Thus, together with 
the standard (DFW) version, we have three distinct versions of the Dirac 
equation in a curved spacetime. The basic quantum mechanics was studied 
for each of those three equations [5] (see Ref. [H] for a summary): 

• the Dirac probability current was defined for any possible field (7^) of 
Dirac matrices. We derived the condition which has to be satisfied by 
(7^) in order that the current conservation be true for any solution of 
the Dirac equation (for DFW, this condition is always satisfied); 

• the precise form of the relevant scalar product was shown to be imposed 
by the axioms of quantum mechanics (it turns out that, for DFW, this 
is the form used by Leclerc [TO]); 

• the Hamiltonian operator was written and the condition of its hermitic- 
ity w.r.t. the relevant scalar product was derived in a general setting 
(for DFW, this condition coincides, in the particular case envisaged by 
Leclerc [10], with that derived by him). 

That foregoing work showed, in particular, that the hermiticity of the Hamil- 
tonian is not stable under all admissible changes of the field of Dirac matrices. 
This implies that there is a non-uniqueness problem for the curved-spacetime 
Dirac equation. 

The aim of the present work JT^, therefore, was to study the (non-) uniqueness 
of the Hamiltonian and energy operators, including the (non-) uniqueness of 
the energy spectrum. The qualitative conclusions are the same for the three 
versions, i.e., the non-uniqueness applies to the alternative equations of Ref. 
[7J, too [TJ. Thus, finding this non-uniqueness was disappointing for us as 
well. In this paper, we will focus on the standard equation (DFW) and the 
non-uniqueness results for it, because these results could be more directly of 
interest to many physicists. 
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2 Three Dirac equations in a curved space- 
time 



The three versions of the gravitational Dirac equation have the same form: 

= -imif), (1) 

where 7 M = 7 M (X) (/i = 0, 3) is a field of 4x4 complex matrices, depending 
on the point X e V in the curved spacetime (V, g^), such that 

7 Y + 7Y = 2ri4, (i,ue{0,...,S} (2) 

[here (g^ v )(X) is the inverse of the matrix (g^^X) and I4 = diag(l, 1, 1, 1)]; 
where ip is a bispinor field for the standard equation (DFW), but is a J^-vector 
field for the two alternative equations, based on the tensor representation of 
the Dirac fields (TRD); 

and where is a covariant derivative, associated with a specific connection. 
For DFW, this is the "spin connection", which depends on the (7**) field 

urn- 

3 Definition of the field of Dirac matrices 

For DFW, one defines 

7"P0 = a»JX) 7« Q , (3) 

with u a = a fi a (X)d[ 1 an orthonormal tetrad field (<9 M is the natural basis 
associated with the coordinates x^) and (7"°) a set of "flat" Dirac matrices, 
i.e., a constant solution of the anticommutation relation (J2J) with (g^) = 
(n Ml/ ) = diag(l, —1, —1, —1). One should be able to use any possible choice 
of (7"°). One should study the influence of both choices: (7"°) and (u a ). 

For TRD, a tetrad field can also be used, though other possibilities also 
exist, e.g. parallel transport [7]. 

In order to be able to use any possible field (7^) of Dirac matrices, we 
must have recourse to the hermitizing matrix A of Bargmann and Pauli [32] . 
This is a 4x4 complex matrix such that 

A* = A, {Afjt = Af /i = 0,...,3, (4) 
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with = M* T = Hermitian conjugate of matrix M. We proved [IT] the 
existence of A, and also that of B: the latter is a positive-definite hermitizing 
matrix for the matrices which are defined by Eq. f JT0|) below. In general, 
A = A(X) is also a field. However, when we use the definition from a tetrad 
field ([3]), we may take [S] 

A = A«, (5) 

where ^ is a (constant) hermitizing matrix for the constant set (7"°). More- 
over, in practice for DFW, one considers only sets (7^) for which 7" is a 
hermitizing matrix. (This includes the usual sets: Dirac's, "chiral", Majo- 
rana's, and any set deduced therefrom by a unitary transformation [T3].) 
Thus, in practice, 

A = 7 fl0 for DFW. (6) 



4 Local similarities 

In a curved spacetime (V, g^u), the Dirac matrices 7^ and the hermitizing 
matrix A are fields: they depend on I 6 V. If one changes from one field 
(7^) to another one (7 M ), also satisfying the anticommutation relation ([2]) 
and with the same metric g^, then the new field obtains by a local similarity 
transformation [12], which is unique up to a complex factor jllj : 

3S = S(X) G GL(4, C) : = S"V(I)5, /i = 0,...,3. (7) 

For the standard Dirac equation (DFW), the similarities are restricted to the 
spin group Spin(l,3), i.e., they result from a change of the tetrad field by 
a (proper) local Lorentz transform L(X) G S0(l,3), depending arbitrarily 
on I 6 V |. The corresponding local similarity transformation S(X) G 
Spin(l, 3) is deduced [8] from L(X) through the spinor representation (defined 
only up to a sign), L i->- S — ±S(L). 

5 The general Dirac Hamiltonian 

Rewriting the Dirac equation $1]) in the "Schrodinger" form: 

^^ = H^, (t = x°), (8) 
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gives the Hamiltonian operator: 



H = ma° -iaWj-iiDo-do), (9) 

with 

a ^7°^ 00 , ^=7°77<7° (7 = 1,2,3). (10) 

One should realize that, in a given spacetime (V, g^v) and with a given field 
of Dirac matrices (7^), this operator still depends on the coordinate system, 
or, more exactly, on the reference frame [S] — the latter being understood 
here as an equivalence class F of local coordinate systems (charts) on the 
spacetime V, modulo the purely spatial transformations 

x 10 = x°, x ' j = f j ((x k )). (11) 

In the present work, we fix the reference frame F and even the chart (the 
latter being, however, an arbitrary one), and we study the dependence of the 
Hamiltonian and energy operators on the field (7^). 



6 Scalar product and equivalent operators 

The Hilbert scalar product is fixed by the following result [Sj: 

Theorem. A necessary condition for the scalar product of time-independent 
wave functions to be time independent and for the Hamiltonian H to be a 
Hermitian operator, is that the scalar product should be 

(iff I <P) = J V'^tV V^9 d 3 x. (12) 

For DFW, we have in practice A = 7" , Eq. ([6]). Then (I12jl is the scalar 
product used by Leclerc [TO] . 

With each of any two possible "coefficient fields": (7^, A) and (7 m ,t4), 
corresponds thus a unique scalar product. [For DFW, we have always A = 
A = Ar, Eq. (J5J, so the coefficient fields reduce to the field of Dirac matri- 
ces 7 M .] These two scalar products are isometrically equivalent through the 
mapping ip H- ip = S~ l ip: one shows easily [1] that 

(f\<f) = I <P). (13) 
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Moreover, H is fully determined by the set of the products (H^ | (p), for 
ip,tp £ T> = Dom(H). The same is true for the Hamiltonian H after the 
similarity ((?]), H being obtained by substituting 7^ for 7^ in Eq. ffTOl) . and 
by using the covariant derivatives = + Y ^ in ([9]), with the new 
connection matrices [see Eq. ffT5"j) below]. It follows [I] that, in order that H 
and H be equivalent operators, it is necessary and sufficient that 

H = S , - 1 HS', (14) 

where S is the relevant similarity in Eq. (T7])- This, of course, is no surprise. 
Of course also, the same condition defines the equivalence of the energy 
operators E and E, E being defined by Eq. (TTTj) below. 



7 Invariance condition of the Hamiltonian un- 
der a local similarity (DFW) 

When does a local similarity S(X), applied to the field of Dirac matrices 
(7^), leave the Hamiltonian operator H [equation ©] invariant? I.e., when 
do we have Eq. (TTJJ? 

The Hamiltonian fl9]) involves the connection matrices, = D^ — d^. For 
DFW, these matrices change after a similarity according to [6]: 

T^S^T^S + S-^d^S). (15) 

A straightforward calculation shows then that, for DFW, we have (Tl4|) iff 
S(X) is time-independent, 

d S = 0. (16) 

But, in the general case, we have 7^ 0, hence any possible field (7^) 
depends on t = x°. Thus, there is no way of finding a class of fields (7^) 
exchanging with OqS = 0. I.e.: The Dirac Hamiltonian is not unique PQ. 
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8 Invariance condition of the energy operator 
(DFW) 

When the Hamiltonian H is not Hermitian, one should use the Hermitian 
part of H : 

E = E+^=B- 1 d {^gB) =i(H + H*), B = Aj°, (17) 

where the Hermitian conjugate operator H* is w.r.t. the unique scalar prod- 
uct (fl"2]) . This Hermitian-symmetrized operator E coincides, in the particular 
case envisaged by Leclerc [10], with what is called the energy operator, which 
is derived from the field Lagrangian. In the general case (thus for TRD and 
for DFW as well, and with any possible (7^) field), the Dirac equation (Op) 
also derives from a Lagrangian [1] . We can show generally that the expected 
value of the energy operator E equals the classical field energy as derived 
from that Lagrangian. 

Again a straightforward calculation gives the invariance condition of E 
(for DFW): 

B(d S)S~ 1 - [B(d S)S~ l ] ] = 2 [B(d S)S~ 1 ] a = 0. (18) 

Only very particular local similarities S(X) do verify (fTBl . Thus, there is 
a serious non-uniqueness problem for DFW (and for the alternative, TRD 
equations as well). Could even the spectrum of E be non-unique? In what 
follows, we investigate this question. 



9 Explicit expression of the energy operator 
(DFW) 

The general expression of the change of E after a local similarity S(X) is 
easily found to be [I]: 

5E = SES" 1 - E = -iB' 1 [B(doS)S~ 1 } a . (19) 

In very general coordinates, the matrix (a M a ) of that tetrad which defines the 
starting (7^) field may be chosen to satisfy a°j = PQ. Then a° = \f(f® 
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from the orthonormality of the tetrad. Let us take for "flat" matrices 7"° 
any usual Dirac matrices, for which A = 7^° [Eq. (jSJ)]. Thus 

B= A 7 = 7 tt Vo7 tt °) = y/j">U, (20) 

whence from ffTUj) : 

5E = -i [(d S)S~ 1 ] a . (21) 

For DFW, S(X) is an arbitrary Spin (1, 3) transformation. Thus, (doS)S~ 1 is 
the generic element of Q, the Lie algebra of Spin(l, 3) - whose the matrices 

s af} = [7 flQ ,7 tt/3 ] (a < P) (22) 

make a basis. Hence 

o"E = -i [uj a ps a ^]\ (23) 

where, depending on the local Lorentz L(X) that defines S(X) = S(L(X)), 
the six coefficients u a /3 = —cup a depend arbitrarily on X G V. With usual 
Dirac matrices verifying ((6]), we have 

(5^)+ = 7 »°s^ T «°, (24) 

whence from (12 3 p and the anticommutation formula: 

3 

5E = -1 u jk s jk . (25) 
j,k=i 

10 The case with the "chiral" Dirac matrices 

If the "flat" Dirac matrices 7"° are the "chiral" ones {see e.g. Schulten (14] . 
Eqs. (10.257) and (10.260)}, we get from ([25]): 

5E = UjkS f = £), JV = -±£* (26) 

i,fc=i ^ ' 

where 6* = (# fc ) with 9 1 = w 2 3 (circular), and where a = (cr^) with a k the Pauli 
matrices. 

Depending on the three real numbers Ujk, 1 < J < k < 3, the matrix iV 
can be any Hermitian matrix 2x2 with zero trace {see e.g. Schulten [14] , 
Eq. (5.226)}. Any such matrix has two eigenvalues \x G R and —fi, and has 
an orthonormal basis of eigenvectors: respectively u G C 2 for /i, and v for 
-/x. 
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11 Non-uniqueness of the energy spectrum 
(DFW) 

A small perturbation of the starting Dirac matrices 7 M is defined by a simi- 
larity close to the identity: 

S(e,X) = I + e(5S)(X) + 0(e 2 ). (27) 

It modifies each eigenvalue of the energy operator E according to the well- 
known formula 

5X = (ip | 6E(e)ip) + 0(£ 2 ), (28) 

where ip is the eigenfunction for the unperturbed state. Inserting <5E (J26]) 
and the expression (TT2"]) of the scalar product into Eq. ()28p . and decomposing 
the bispinor into two two-components spinors: ip = (0, x), we get: 

5\ = J ip j SEip ^-gg 00 d 3 x = J (0 f iV0 + x ] Nx) dV. (29) 

Let us fix fj, > and t. For any spatial position x = (x J ) in our fixed chart, 
let us select the 2x2 Hermitian matrix with zero trace, N = N(x), in such 
a way that 0(x) be the eigenvector of N(x) for the eigenvalue fi. (The 3- 
vector x is the coordinate expression of the position x in the space manifold 
M associated with the arbitrary given reference frame F |S].) That is: 

$N(f> = n (t) ] (f). (30) 

Since N has the eigenvalues \x and — fi, we have then necessarily 

X ] N X > -fi X ] X- (31) 

In Eq. (13 ip . the inequality becomes an equality only if x( x ) is an eigenvector 
of N(x) with eigenvalue —fi. In that case, 0(x) is orthogonal to x( x ) i n C 2 , 
x(x)±0(x). Returning to Eq. ( )29|) . we see that 5\ > unless if i) x( x )-L0(x) 
for almost every spatial position x, and ii) J <fr(j) dV = J x^X dV. This could 
occur only in extremely particular situations. In fact, i) alone implies pQ that 
the probability current is light-like for almost every x, and this is impossible 
if m > [T5]. This proves the non-uniqueness of the DFW energy spectrum, 
at least for a massive particle. 
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12 Conclusion 



Three distinct gravitational Dirac equations were envisaged: the standard 
equation ( "Dirac- Fock-Weyl" or DFW), plus two alternative equations [7j 
based on the tensor representation of the Dirac field (TRD). Using the her- 
mitizing matrix A, they were studied together, to a large extent [Eld]. Non- 
uniqueness results have been proved for each of the three versions pQ. In this 
paper, we summarized the latter work, but focussed on the results for DFW. 

The Hamiltonian operator H is not unique: it depends on the admissible 
choice of the field of Dirac matrices. The same is true for the energy operator 
E. This is true for DFW and for TRD. 

The spectrum of E is itself non-unique. All of these results apply already 
to a flat spacetime in a non-inertial frame. 

A Appendix: The classical energy and its 
frame dependence 

One might ask if the non-uniqueness of the energy spectrum for a spin | par- 
ticle in a curved spacetime, reported here, might have something to do with 
the well-known fact that, in general relativity (GR), there is no covariant 
concept of local energy — already in the classical context (as opposed to the 
quantum context envisaged in this paper). The latter fact may be formulated 
by saying that there is no covariant conservation law for the (total) energy- 
momentum — the energy-momentum of the gravitational field giving rise only 
to a pseudo-tensor. As is well known, the usual "conservation law" for the 
energy-momentum tensor of matter and non-gravitational fields "does not 
express the conservation of anything" [16] . However, for a classical test par- 
ticle, in any arbitrary reference frame F, there is a well-defined Hamiltonian 
energy [TJ [T71 [TH], though it depends on F and, in general, on the time: 

• Geodesic motion in the Lorentzian manifold (V, g^ u ) derives from the 
("super-") Hamiltonian over the 8-dimensional phase space associated 
with V: 

H[(P»), CO] = \r{{x P ))P»P» (c = I)- (32) 

This is because, in a (pseudo-)Riemannian manifold, geodesic motion is 
identical to free Hamiltonian motion, i.e., H = T with T[(p M ), (x v )] = 
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^Q^PfiPu the "kinetic energy": Sect. 45 D in Arnold [TJJ]. Note that H 
does not depend on the (proper) time r. 

• It follows from that "time" -independence that, in any given coordi- 
nate system, we deduce a "normal" Hamiltonian H over the usual, 
6-dimensional phase space, by dimensional reduction (Sect. 45 B in 
Ref. [19]): 

H[( Pj ),(x^,t = x°]= Po (33) 

is extracted from 

gTp&v ~ m 2 = 0. (34) 

{With the ( — h ++) signature, we have H = — po instead, and the 
last equation writes g^p^Pu + m 2 = [E].} One may note that E = 
H = po is invariant under purely spatial coordinate changes (TTTT) . thus 
it depends only on the reference frame — which is arbitrary, but fixed in 
the present study. 

Therefore, in a fixed reference frame, a classical test particle in a given exter- 
nal gravitational field has unique energy, in spite of the absence of a covariant 
concept of local energy in GR, where gravitational dynamics are also included 
in the model. Similarly, the answer to the question posed at the beginning of 
this Appendix is no. The non- uniqueness of the energy spectrum for a spin- 
half particle in a given external curved spacetime, as found in this study, has 
nothing to do with the absence of a covariant concept of local energy in GR. 
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